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Precise understanding of non-linear evolution of cosmological perturbations during inflation is 
necessary for the correct interpretation of measurements of non-Gaussian correlations in the cosmic 
microwave background and the large-scale structure of the universe. The "8N formalism" is a 
popular and powerful technique for computing non-linear evolution of cosmological perturbations 
on large scales. In particular, it enables us to compute the curvature perturbation, £ , on large scales 
without actually solving perturbed field equations. However, people often wonder why this is the 
case. In order for this approach to be valid, the perturbed Hamiltonian constraint and matter-field 
equations on large scales must, with a suitable choice of coordinates, take on the same forms as 
the corresponding unperturbed equations. We find that this is possible when (1) the unperturbed 
metric is given by a homogeneous and isotropic Friedmann-Lemaitre-Robertson- Walker metric; and 
(2) on large scales and with a suitable choice of coordinates, one can ignore the shift vector (goi) 
as well as time-dependence of tensor perturbations to gij/a 2 (t) of the perturbed metric. While the 
first condition has to be assumed a priori, the second condition can be met when (3) the anisotropic 
stress becomes negligible on large scales. However, in order to explicitly show that the second 
condition follows from the third condition, one has to use gravitational field equations, and thus 
this statement may depend on the details of theory of gravitation. Finally, as the SN formalism 
uses only the Hamiltonian constraint and matter-field equations, it does not a priori respect the 
momentum constraint. We show that the violation of the momentum constraint only yields a 
decaying mode solution for (, and the violation vanishes when the slow-roll conditions are satisfied. 



I. INTRODUCTION 

Given the success of cosmological linear perturbation 
theory, the focus has shifted to non-linear evolution of 
cosmological perturbations. As the magnitude of the 
primordial curvature perturbation is of order 10 -5 , any 
non-linearities are expected to be small; however, such 
non-linearities can be measured using non-Gaussian cor- 
relations of cosmological perturbations (such as tempera- 
ture and polarization anisotropy of the cosmic microwave 
background [l[ and density fluctuations in the large-scale 
structure of the universe @] ) ■ For this reason, precise un- 
derstanding of the non-linear evolution of cosmological 
perturbations is of great interest in cosmology. 

The so-called U 5N formalism" 0-0] is a popular tech- 
nique for computing non-linear evolution of cosmological 
perturbations on large scales. Here, by "large scales," 
we mean the scales greater than the Hubble horizon, in 
a sense that the comoving wavenumber of perturbations, 
k, is much less than the reciprocal of the comoving Hub- 
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ble length, i.e., k <C aH. In particular, it enables us to 
compute the curvature perturbation, £, without actually 
solving the perturbed field equations. In this paper, we 
show why this is the case by re-deriving the SN formal- 
ism using the gradient expansion method as applied to 
Einstein's field equations and scalar-field equations in the 
flat gauge. The usual derivation of the SN formalism is 
based on the so-called "separate universe" approach 
which assumes the existence of a locally homogeneous 
(but not necessarily isotropic region smoothed over 
some large length scale. We provide a support for this as- 
sumption by considering a global region including many 
such smoothed local regions, and show that they behave 
as if they were locally homogeneous regions which evolve 
independently from each other. In so doing, we point out 
a subtlety regarding the momentum constraint, which is 
not a priori respected by the SN formalism. 

The organization of this paper is as follows. In Sec. II, 
we describe our basic setup including the metric, gauge, 
and scalar-field Lagrangian. In Sec. Ill, we review the 
gradient expansion method, which constitutes the basis 
for the SN formalism. In Sec. IV, we re-derive the SN 
formalism. In Sec. V, we give the sufficient conditions 
for the validity of the SN formalism and conclude. 
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II. BASIC SETUP 



A. Metric 



We write the spacctime metric in the Arnowitt-Deser- 
Misner (ADM) form, which is the standard (3+1)- 
decomposition of the metric [l(| , 

ds 2 = -a 2 dt 2 + (dx l + /3 l dt) (dx 3 + [3 3 dt) , (1) 

where 7^ is decomposed as follows, 



7^aV* (e h ) i 



(2) 



Here, a is the scale factor which depends only on time, 
and ip is the scalar perturbation to the spatial curvature. 
The traceless tensor, hij, is further decomposed as 



hij = diCj + djd - -Sijd k C k + h 



(T) 
ij ) 



(3) 



where Cj contains both scalar and vector perturbations, 

(T) 

whereas represents tensor perturbations. 

We decompose the extrinsic curvature, Kij, into a 
trace part, K, and a traceless part, Aij, atQ 



Hi 



K + a 2 e 2 ^A t 



(4) 



Einstein's field equations written in terms of these vari- 
ables are summarized in Appendix A. 



B. Flat gauge 

In this paper, we shall fix the gauge completely (i.e., 
leaving no gauge degree of freedom) by imposing the fol- 
lowing gauge-fixing condition^: 



•0 = d = 0. 



(5) 



Therefore, the spatial metric is described only by the 



scale factor and tensor perturbations as 7^ 



This gauge was also used by [11| (see his Eq. (3.2)). 

We shall call this gauge the "flat gauge" throughout 
this paper. Note that a flat-gauge condition in the liter- 
ature sometimes does not include Cj = 0. In such a case 
a residual gauge degree of freedom would remain. In the 
flat gauge, the metric is given by 



ds 2 



a 2 dt 2 + a 2 {t)[e h ']y (dx % + p l dt) (dx 1 + p'dt) . 

(6) 



1 In Salopek and Bond [f| , a trace- free part of the extrinsic curva- 



ture is denoted as Ka . In our notation, wc have Ki 



Sometimes another gauge condition, ip = and ft = 0, is im- 
posed; however, the degree of freedom of this gauge is not com- 
pletely fixed, as this gauge condition only gives the relation for 
the time derivatives of the shifts in temporal and spatial coordi- 
nates. 



The variables in this metric such as a, pi, and con- 
tain non-linear perturbations. However, we shall assume 
that the unperturbed metric is still given by a homo- 
geneous and isotropic Friedmann-Lemaitre-Robertson- 
Walker metric: 



ds 2 



-dt 2 + a 2 (t)S ij dx i dx j . 



(7) 



Therefore, our argument below does not hold if the un- 
perturbed metric is not given by Eq. (jTj) . 



C. Scalar-field Lagrangian 



We shall consider a universe filled with scalar fields: 



C = --G I jd^ I d fi i P J -V. 



(8) 



The capital Latin indices (/, J, etc.) denote scalar-field 
components running from 1 to n where n is the number 
of scalar fields. Here, Gjj is the metric tensor for scalar- 
field space. For simplicity we shall take the canonical 
kinetic term, G/j = Sjj, for the moment. We then argue 
later (in Sec. lIII Cj) that the results are also valid for non- 
canonical kinetic terms in the first order of the gradient 
expansion. 

With this Lagrangian, the stress-energy tensor and the 
field equation of scalar fields are given by 

Tfi V = Gud lt <p I 8 v <p J + \ -Gu^ff x <p I d a <p J - V^j , 

(9) 

ygg^dutf) -Vj = 0, (10) 
where Vi = dV/dip 1 . 



III. GRADIENT EXPANSION METHOD 

A. Ordering in the gradient expansion 

Since we are interested in non-linear perturbations on 
super-horizon scales, we shall expand field equations in 
the number of spatial derivatives: this is called the gra- 
dient expansion method 0,0- In this method, the ratio 
of the comoving wavenumber and the comoving Hubble 
scale. 



k 



(11) 



is taken to be a small parameter. 

Before we proceed, let us emphasize that we assume 
the validity of perturbative expansion: namely, while we 
shall deal with non-linear perturbations, we assume that 
the (i + l)th-order perturbations are smaller than the 
ith-ordcr perturbations. This means that we have two 
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smallncss parameters: one is the number of derivatives, 
e, and the other is a smallness parameter of perturbation 
theory, S, which corresponds to ip, /3 ! , ip 1 — (p 1 , etc. 

These two parameters should satisfy the following con- 
dition: 

6<e<l. (12) 

This is because, if we take the smoothing length to be 
infinitely large, i.e., e — ¥ 0, then the perturbation must 
vanish, i.e., 5 — > 0. Then, the metric must approach the 
unperturbed metric given by Eq. ([7]) as we take e — > 0. In 
other words, the amplitude of the perturbations should 
be limited by the smoothing length we take. 

We now estimate the ordering of perturbation variables 
in terms of the gradient expansion. First, we demand 
that all physical quantities do not vanish in the lowest 
order of the gradient expansion: 

a - 1 = Jij - 8 l3 = ip 1 - ip 1 = O(e , S), 

F = o( < r\S). (13) 

We do not include ijj here because we work in the flat 
gauge. We have defined jij = jij /[a 2 {t)e 2 ^], which is 

equal to [e h ]ij in the flat gauge. 

One may wonder why we chose to start with (3 l = 
0(e^ 1 ), which seems to diverge in the limit of e — > 0. 
However, this is not true. As noted earlier, the existence 
of the perturbation (S > 0) guarantees e > and we 
always have (3 l = 0{e-\S) < 1 from Eq. (12]); thus, 
there is no divergence in the metric. In fact, we recover 
the standard Friedmann equation in the lowest order ap- 
proximation. Furthermore, at the end of Sec. IIVD1 we 
show that consistency between the Hamiltonian and mo- 
mentum constraint equations demands /3 l = C(e _1 ). 

Note that the shift vector comes with a spatial deriva- 
tive, di, in Einstein's field equations and scalar-field equa- 
tions. As dif3 l — 0(e°), the spatial derivatives are kept 
in Einstein's field equations and scalar-field equations for 
a — 1 = jij — Sij — ip 1 — (p 1 = 0(e°). In other words, 
as we keep spatial derivatives in our approach, we are 
considering some global region in which there are many 
smoothed local regions. Therefore, we do not a priori 
demand that these local regions evolve independently of 
each other, contrary to what is always demanded by a 
separate universe approach. Specifically, for a separate 
universe approach, j3 l = 0(e) is always assumed a priori. 

Similarly, when we decompose the quantities Pi and 
Ci into scalar and vector components as Pi = diP^ s > + 
P^ and Ci = diC^ + C± , respectively, the scalar 
components arc of order e -2 : (3^ = C' 5 -' = 0(e~ 2 ). 

In order to see how Eq. (| 13[) can be relaxed, we now 

investigate the nature of solutions for Pi and h^J . 

At the first order in perturbation variables and the low- 
est order in the gradient expansion, the evolution equa- 
tion for Aij is given by (see Eq. (|A10[0 : 

ly+3ifly = 0(e,rj 2 ), (14) 



where H is the Hubble expansion rate, H = a/a. Here, 
we have ignored the anisotropic stress term on the right 
hand side of Eq. (|A10|) . as it is of the second order in the 
gradient expansion for scalar fields. This is a stronger- 
than-necessary condition: Eq. (fT4|) is still valid if the 
anisotropic stress of matter fields is of the first order in 
the gradient expansion. 

It follows from Eq. (fl"4"]) that the traceless part of the 
extrinsic curvature, Aij, has a decaying solution, Aij cx 
1 jo? [|| On the other hand, the evolution equation for 
jij with d = yields (see Eq. ([A"8|) 1: 

h$p - -2Aij + ^ (dipj + djfo - ^5ij6 kl d k pij . (15) 

As the scalar, vector, and tensor modes arc independent 
in linear theory, the equations for the shift vector and 
tensor perturbations arc given by 

+3Hh ( i f ) =<D{e,5 2 ), 
/3 i + 3H(3 i =0(e a ,6 2 ). (16) 

Therefore, h[J^ and /3 4 = /J; /a 2 also have decaying so- 
lutions scaling as a~ 3 . This is a consequence of the fact 
that the unperturbed metric (Eq. (0 ) is given by a homo- 
geneous and isotropic Friedmann-Lemaitre-Robertson- 
Walker metric. In other words, this result may not hold 
for anisotropic models such as Bianchi-type metrics. 

At the second order in perturbation variables, as the 
source terms in Eq. ()A8j) and Eq. (|A10|) are decaying, 

(T) 

the second-order equations for and fii have ap- 

proximately the same forms as the first-order equations 
(Eq. (fT6)) and Eq. (fT4]) ), and thus their solutions must 
also be decaying as a~ 3 . Similarly, n-th order solutions 
for n > 3 are also decaying. 

These properties allow us to safely ignore, in the low- 
est order of the gradient expansion and the n-th order 
of perturbation theory, the traceless part of the extrin- 
sic curvature A^, the shift vector p\ as well as a time 

derivative of tensor perturbations, , after the decay- 
ing solutions become sufficiently small. This means that 
these quantities must be higher order in the gradient ex- 
pansion than naively assumed in Eq. (|13|) : Aij = 0(e), 
=O(e ), and h\j ] = 0(e). 

Now, it turns out that the above argument also applies 
to the next order of the gradient expansion. At the next 
order in the gradient expansion: 

a-l = li j-8ij = i P I -(p 1 = 0(e), 

P l = 0(e°), (17) 



3 As we start with /J 1 = 0(e 1 ), we need to linearize Eq. d AlOl l 
to obtain Eq. (1141 . showing Aij cx 1/a 3 . On the other hand, 
assuming fi l = 0(e), Hamazaki derives AijA % i oc 1/a 6 without 
using perturbation theory (see Eq. (2.54) of [T3l ). 



4 



one can show that, for scalar fields whose anisotropic 
stress is of the second order in the gradient expansion, 
the equations take on the same form as Eq. (fT4|) and 

Eq. or 



Aij 

m _ 



3Hh^=0(e 2 ,S 2 ), 



ft + 3#ft = 0(e,<5 2 ). 



(18) 
(19) 
(20) 



Applying the same argument as above, one finds that 
Aij , h\j , as well as ft decay for the n-th order in per- 
turbation theory. We thus find: A tJ = 0(e 2 ), ft = 0(e), 

and h^p = 0(e 2 ). However, this argument cannot be ex- 
tended to the second order of the gradient expansion, as 
Eq. (|18p is valid only when the anisotropic stress term is 
unimportant. As one can no longer ignore the anisotropic 
stress of scalar fields at the second order in the gradient 
expansion, Eq. (|18|) is no longer valid in that order. 
Therefore, Eq. ([TUf should be revised as 



a- 1 
ft = 0(e), h 



9 

(T) 



<^ = 0(e°), 
0(e 2 ), 



(21) 



where we have dropped 5 in 0(. . . ), as the above estima- 
tion is valid for all orders of perturbation theory. Note 
that this result is valid only in the fiat gauge given by 
Eq. ([5]). In particular, the condition d = was needed 
to estimate the gradient-expansion-order of ft and . 

These results might depend on the details of theory 
of gravitation, as we have used Einstein's field equations 

' (T) 

to obtain solutions of ft and h\- . Furthermore, as a 

' (T) 

perturbative expansion is used in estimating ft and h\- , 
the validity of a perturbative description of the metric 
with the unperturbed metric given by Eq. ([7]) has been 
assumed in the above argument. 



B. Comparison with previous work 

How does Eq. (pH]) compare with the previous work? 
Our starting point, Eq. (|13[) . is different from the assump- 
tion made in Lyth, Malik and Sasaki [g| (also see jl3l|). 
They assume that there exist an approximate set of co- 
ordinates with which the metric of any local region can 
be written as a Friedmann-Lemaitre-Robertson- Walker 
metric. This implies that the shift vector, ft, vanishes 
and the quantity, 7^, is time- independent in the limit 
of e ->• 0: ft = 0(e) and %j = 0(e). We do not make 
this assumption a priori, and thus our argument is more 
general than that given in They then show that, us- 
ing Einstein's field equations and ignoring the anisotropic 
stress term, 7^ decays in the first order of the gradient 
expansion, concluding 7^ = 0(e 2 ). 

In [l4j . Weinberg uses a broken symmetry argument 
to show ft = 0(a~ 2 ) and 7^ = 0(a~ 2 ) for generally 



covariant theories and with a suitable choice of coordi- 
nates, assuming that the unperturbed metric is given by 
Eq. ([7]) and the anisotropic stress term is negligible. He 
then shows that, for Einstein's field equations in a co- 
ordinate system in which ft = and a certain combi- 
nation of matter perturbations vanishes, this solution is 
an attractor. By identifying 0(a~ 2 ) with 0(e 2 ) (because 
each spatial derivative must come with 1/a), his argu- 
ment yields ft = 0(e 2 ) and 7^ = 0(e 2 ). 

Therefore, our finding agrees with the previous work: 
if the unperturbed metric is given by Eq. ([7|), the 
anisotropic stress term is negligible on large scales, and 
field equations are given by Einstein's field equations, 
then ft = 0(e) and % = 0(e 2 ). Note that Weinberg's 
estimate for the order of ft is higher by e; however, he 
does not show ft = 0(e 2 ) explicitly because he chooses 
coordinates in which ft = 0. 



C. Gradient expansion of Einstein's field equations 
and scalar-field equations in the flat gauge 

Now, we apply the gradient expansion to Einstein's 
field equations and scalar-field equations. We shall work 
with the flat gauge given by Eq. ([5]), which gives the 
gradient-expansion order of perturbation variables given 
in Eq. (EJ). 

We shall choose the number of e-folds, N = J. Halt' , 
as our time coordinates in the flat gauge. The Hamilto- 
nian constraint (Eq. (|A5[) ) and the scalar- field equation 
(Eq. (jTDJ) ) in both the lowest order and the next order of 
the gradient expansion are given by 



3H 2 M 2 



P, 



(22) 



Hd N (Stplr) + 3H 2 ^ N + Vi = 0, (23) 

where H = H/a is related to a trace of the extrinsic 
curvature as H = —K/3, and the subscript N denotes a 
partial derivative with respect to N. The energy density, 

p, is given by p — ^-Gijip^tf^ + V. Then, H is given 

by m 



H 2 



2V 



6Af 2 - Gu^ nV j n 



(24) 



coupled with 



All we need to do is to solve Eq. 
Eq. (El 

On the other hand, the unperturbed equations are 



3ff 2 M 2 = p, 



(25) 



Hd N (H^ N ) +3H 2 <p I N + V I (£) = 0, (26) 

where p and (p 1 are the unperturbed energy density and 
scalar fields, respectively. Apparently, the perturbative 
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equations (Eqs. (|22|) and (|23[0 coincide exactly with the 
unperturbed equations (Eqs. (|23|) and (f2l))) ). This re- 
sult shows that each region smoothed by a super-horizon 
scale, e 1, in the universe evolves independently and 
behaves like an unperturbed universe, providing a sup- 
port for the assumption made by a separate universe ap- 
proach. 

These results might depend on the details of theory of 
gravitation. While the correspondence between the per- 
turbed and unperturbed equations for other theories of 
gravitation is an interesting problem, in this paper we 
shall focus on Einstein's General Relativity. However, 
these results should not depend on the form of the La- 
grangian of scalar fields. This is because the anisotropic 
stress (i.e., a traceless part of the stress-energy tensor) for 
scalar fields with arbitrary Lagrangian necessarily comes 
with two spatial derivatives, and thus it must be C(e 2 ). 
As a result, A{j has a decaying solution and a Fricdmann- 
Lcmaitrc-Robcrtson- Walker universe will be restored on 
large scales. 

One can generalize the above results for the canon- 
ical case to non-canonical Lagrangians given by C = 
P{X IJ ,tp K ), where X IJ = -g^d^d^tf 1 . The Hamil- 
tonian constraint is still given by Eq. (|22|) while the scalar 
field equation is given by 



Hd 



A' 



Hd 



n<Pn 



p IJ +H 2 v J N d N p I j+'m 2 p I jv J N +^- = o, 



where P u = 8P/dX 1J , X 
ergy density is defined as p 



I J 



(27) 

= H 2 ip I N ip J N: and the en- 
2X IJ P U -P. As P and 
, and H, one can write H 
if an explicit form of the 



Pu are functions of X IJ , 
as a function of X IJ and 
Lagrangian P is specified. 

What are the implications of these results? As the 
equations take on the same forms, the functional forms 
of the solutions for the perturbed equations and those for 
the unperturbed equations must be the same. Therefore, 
the perturbed solutions, (p 1 , are given by the unperturbed 
solutions, (p , with perturbed initial conditions computed 
in the flat gauge, <pl(x) and ip^^x): 



y I (N,x) = <p I (N,ipi{x),^{x)). 



(28) 



This is the fundamental result of the gradient expansion 
as applied to Einstein's field equations and scalar-field 
equations in the flat gauge. Here, the subscript * indi- 
cates that the quantity is evaluated at some initial time, 
where all the relevant fields are sufficiently outside their 
sound horizon, i.e., k <C a(t*)H(t*) / 'cf, where is the 
speed of sound of propagation of an 7th scalar field per- 
turbation. 

In order to simplify our notations, from now on we shall 
use the lower-case alphabet indices, such as a, b, c . . . , 
to denote the numbers of scalar fields and their time- 
derivatives: 

with a running from 1 to 2n. With this notation, the solu- 
tion (Eq. ([25])) is expressed as ip a (N, x) = <f a (N, fl{x)). 



Similarly, we can write the perturbed energy density of 
multi-scalar fields using the unperturbed energy density 
solution: p(N,x) = p(N,ip"(x)). 



IV. THE 5N FORMALISM 

We now need to relate perturbed initial scalar fields 
and their derivatives, (p%(x), to the observables. In cos- 
mology, it is now customary to express the observables 
such as temperature and polarization anisotropies and 
the large-scale distribution of galaxies in terms of a cur- 
vature perturbation in the "uniform-density gauge," de- 
noted as C- 

The so-called SN formalism 0-0] achieves this by re- 
alizing that £ is equal to a perturbation to the number 
of e-folds, N, arising from perturbed initial scalar fields, 
ip%(x), computed in the flat gauge. 



A. Conservation of £ outside the horizon 

We define the uniform-density gauge as Sp = and 
C{ = 0. (Once again, a uniform-density gauge in the 
literature sometimes does not include Cj = 0. In such a 
case the gauge is not completely fixed.) Let us denote a 
value of tp in the uniform-density gauge as V , |5jo=Ci=cH 
and write £ as 



C = 1p\sp=Cf. 



(30) 



This quantity is useful for extracting information 
about the physics of inflation, as it is conserved out- 
side the horizon, provided that the adiabatic condition, 
p = p[p], is satisfied [f| Q. This is easily seen from the 
energy conservation equation in the lowest order as well 
as in the next order gradient expansion with the gauge 
condition C; = 0: 



p + 3(H + i>\ Cz=0 )(p + p) = 0. 



(31) 



The perturbation to this equation in the uniform-density 
gauge yields 



ip\sp=d=o = 0, 



(32) 



and thus C = V'|5p=c i =o becomes a constant, provided 
that the adiabatic condition is satisfied. 



Incidentally, as the gauge is completely fixed for ip\sp=O =0i 
there is no ambiguity with respect to the residual gauge degree of 
freedom. Moreover, one can always write perturbation variables 
(such as after gauge-fixing as a combination of perturbation 
variables before gauge-fixing (such as d and Sp) such that tp is 
explicitly gauge-invariant. For example, we have, at the linear 
order, 

,, , d,C n 8p 

C = V\Sp=Ci=Q ~ " 



(29) 



3 pjv 

where the right-hand-side of Eq. I I29II is the well-known form for 
a gauge-invariant curvature perturbation in the linear order |16|| . 
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Alternatively, Eq. pip may be integrated with respect 
to t without imposing Sp = 0: 



( = 4> + 



dp 



p 3(P + P\p\) 



= const. 



(33) 



One may then identify this quantity, £, as a general- 
ization of C when dp = is not imposed; however, £ 
is not gauge-invariant and does not coincide with £ = 
ip\sp=d=o, unless the adiabatic condition is satisfied. 

What about scalar fields? As scalar fields do not satisfy 
the adiabatic condition in general, £ is not conserved in 
a universe filled with scalar fields. However, as shown by 
, £ is generally conserved outside the horizon when in- 
flation was driven by a single scalar field. More precisely, 
£ is conserved outside the horizon in a universe domi- 
nated by a single scalar field, provided that the slow-roll 
conditions are satisfied, or that we completely neglect a 
decaying mode solution without imposing the slow-roll 
conditions. This implies that the slow-roll conditions 
correspond effectively to the adiabatic condition and the 
neglect of a decaying mode solution for a single scalar 
field. 



B. Relation between £ and the difference in the 
number of e-folds 

The relation between the curvature perturbation and 
the number of e-folds is given by the gauge transforma- 
tion of the spatial metric, Ty. Under a gauge transfor- 
mation given by t -> T = t + ST and x i -> X 1 = x l + 
the metric transforms as gij(t,x) — > gij(T,X). Let us 
write the 3-metric in the original coordinates in terms of 
the 3-metric in the new coordinates: 



a 2 (T,X) 
lu{T,X) 



dST dST 



dx l 




dX k 


dST 


dx l 


dxi 


dX k 


dX l 



dx l dxi 



MT,X) 



dX k dST 
dxi dx % 

(34) 



We shall always impose Ci = 0, which completely fixes 
the spatial gauge degree of freedom, and thus we can set 
= without loss of generality. 

Let us examine each term in terms of the gradient ex- 
pansion order. The first term is of order 0(e ). As shown 
in Sec. 1111 A[ when Cj = 0, the shift vector is of 0(e); 
thus, the second and third terms are of order C(e 2 ). This 
means that, up to C(e 2 ), the 3-metric transforms as 



7«|c,=o(*> as) = %\c t =o{T, x) + 0{e 2 ). 



(35) 



Recalling 7^ = a 2 (t)e 2 ^(e h )ij and taking the deter- 
minant and logarithm of the both sides of Eq. (f35|) . we 
find 



il>\c i =o(t,x) = il>\c i =o(T,x)+1n 



a(T) 
a(t) 



(36) 



Thus, tp\c i= o approximately transforms as a scalar quan- 
tity having ln(a) as the unperturbed value. It follows 
from Eq. (|36| that the gauge transformation of ip\d=o 
from the flat gauge (in which -0 = 0) into the uniform- 
density gauge is given by 



ip\6p=Ci=o(T,x) = In 



a(t) 



(37) 



AT) J' 

where T denotes time coordinates in the uniform-density 
gauge. 

On the other hand, when we go from the flat gauge to 
the uniform-density gauge, the number of e-folds, N = 
J. Hdt' ', transforms as N — > N, where 



N 



7 H(t')dt' = lnf^T 
. \o(t.) 



(38) 



Here, is an arbitrary initial time. 

Comparing Eq. ([37jl to Eq. (|38|) . one finds 



ip 1 5 P =Ci=o (T, x) = In 



a(t) 



N — N = 5N. 



In 



a(T) 
a(t*) 



(39) 



Therefore, C = 4'\8p=Ci=o is equal to the difference be- 
tween the number of e-folds computed in the flat gauge 
and that computed in the uniform-density gauge. The 
remaining task is to relate SN to perturbed initial scalar 
fields in the flat gauge. 



C. Relation between SN and perturbed initial 
scalar fields in the flat gauge 

The most important result that came from the gra- 
dient expansion of Einstein's field equations and scalar- 
field equations in the flat gauge is that perturbed quanti- 
ties can be calculated using their unperturbed solutions 
with perturbed initial scalar-field values and their time 
derivatives computed in the flat gauge. Therefore, a 
perturbed energy density in the flat gauge is given by 
p(N,x) = p(N,(fil(x)). Here, we choose the number of 
e-folds as time coordinates. 

On the other hand, by definition the energy den- 
sity in the uniform-density gauge (whose time coordi- 
nates are denoted as N) is equal to the unperturbed 
density. Namely, when we go from the fiat gauge to 
the uniform density gauge by changing the time coor- 
dinates as N —> N = N + SN, the density transforms 
as p(N,x) — > p(N,x) = p(N). Here, Ci = is satisfied 
in both gauges, and thus there is no ambiguity with re- 
spect to the spatial gauge degree of freedom. Now, as 
the energy density is a four scalar, 



p(N,x) = p(N,x) = p(N), 



(40) 



which gives p(N,(p^(x)) = p(N). Inverting this result 
yields 



N = N(p,<ft(x)), 



(41) 
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where the functional form of N is the same as that of the 
unperturbed number of e- folds. That the unperturbed 
density p (not t or N) is used as the time coordinates 
here ensures that the final time slice coincides with the 
uniform density hypersurface. 

With these results, we can finally calculate f : 

( = N-N 



(42) 



where S(p"(x) = ip^{x) — (p% denotes perturbations to 
initial scalar fields computed in the flat gauge, and N a 
and N a b are defined as 



N„ 



N„ 



d 2 N[p,& 
d(p%d(pl 



(43) 



This is the 8N formalism, which enables us to relate £ to 
the initial scalar field perturbations (i.e., scalar field per- 
turbations at the initial time) computed in the flat gauge, 
once we know derivatives of the number of e-folds with 
respect to the initial values of the unperturbed scalar 
fields, ipl, and their derivatives, </?jv*- 



D. Momentum constraint 

Perhaps a striking thing about the SN formalism is 
that we only had to use the Hamiltonian constraint 
(Eq. p2"])) and the scalar- field equation (Eq. ([121)) hi the 
gradient expansion. But, should not we also impose the 
momentum constraint for consistent calculations? 

As the momentum constraint comes with a spatial 
derivative, we need to consider the momentum con- 
straint in 0(e 2 ) in order to derive the correct relationship 
between physical quantities up to 0(e): 

diH = -^Gjj^d^ + 0(e 3 ), (44) 



where we have used the fact that Aij = C(e 2 ). 

On the other hand, the Hamiltonian constraint 
(Eq. ([22]) ) may be differentiated by di to give 



diH = -JLfGut&dup* + B u (45) 



where 



H 3 

Bi = ^F G u (fNdiip J N - fif^f 7 ) . (46) 

Here, we have used the equation of motion for scalar fields 
given by Eq. (|23|) , as well as the evolution equation for K 

given by Eq. (|A9p . which yields Hm = —^Ti^^'IJ'-Pn'Pn 

p 

in the gradient expansion in the flat gauge. 



Comparing Eq. (J45j) and Eq. |44]) , we find that the SN 
formalism, which does not use the momentum constraint 
but uses only the Hamiltonian constraint, can violate the 
momentum constraint by an amount Bi. 

How important is B{1 In order to investigate the be- 
havior of Bi, let us take a spatial derivative of the equa- 
tion of motion for scalar fields: 



H d 



Vu 



5 AC dN 



+ H^{B lN + 3Bi) + 2(p I NN HB i = 0. 

(47) 

By contracting this equation with ipj^, one finds 



d N B, + 3B t = 0, 

alB„ 
— >• Hi = ^ — ■ 



(48) 



Therefore, Bi has only a decaying solutionis This is a 
good news for the SN formalism: while it does not a pri- 
ori respect the momentum constraint, a violation of the 
momentum constraint rapidly decays away by inflation. 

That Bi is a decaying mode may be traced back to 
the fact that the traceless part of the extrinsic curvature, 
Aij , is a decaying mode in the gradient expansion in the 
flat gauge. In other words, it is a consequence of the 
universe behaving like a Friedmann-Lemaitre- Robertson- 
Walker universe on super-horizon scales, which is guar- 
anteed by Eq. ([211 . 

One may wonder why consistency between the 
momentum constraint and the Hamiltonian con- 
straint gives a relation only among scalar fields: 
Gij (ifj^diipff — (p I NN di(p J ) = 0, rather than a relation 
between the metric variables and scalar fields. This is 
because we have ignored the decaying solutions of (3 l and 

' (T) 

h\j ■ To see this, let us work at the first order in per- 
turbations, and bring f3 l back into Einstein's field equa- 
tions. We find that consistency between the momentum 
constraint and the Hamiltonian constraint gives 



Bi 



M 2 H 2 



(49) 



Indeed, consistency gives a relation between scalar fields 
(contained in Bi) and a metric var iable (/^jj This 
equation also indicates that, when we keep the decay- 
ing quantities of order a -3 , the gradient-expansion order 
of ft is indeed 0(e _1 ) (see Eq. CE3|)), as Eq. dMJ) gives 
dip = 0(1). 



5 In linear theory, this quantity is equal to 
given by Eq. (2.26) of |15| . as well as to 
by Eq. (5.23) of [Jj]. They show W oc a 
theory. 

6 We thank M. Sasaki for clarifying this point 



a~ 3 diW where W is 
—dif where / is given 
3 and / oc a - 3 in linear 



8 



Now is the time to answer the following question: what 
if we demand f3 l = C(e )? In this case Eq. (|4"9"|) gives 
Bi = 0(e 2 ), or 



Gij 



N 



^ NN 5 V J ) = 0(e), 



(50) 



where Sip 1 = Lp 1 — (p 1 is the perturbation of scalar fields 
in linear theory. This result indicates that we are not al- 
lowed to have a configuration of scalar fields which yields 
Gij (tpN^N ~ <P I NN<> i P J ) = ^( e °)- On the other hand, 
if we start with /3 1 = 0(e~ 1 ), then we can show that 
Gij (0 n S<Pn — 'f I NN dip J ) becomes negligible as it is a 
decaying mode. Another way of saying this is that, if we 
start with j3 l = O(e ), then Gij (ip^Sipjf - ip^j^Sip' 1 ) = 
gives only one solution for Stp 1 , and another solution, 
which corresponds to a decaying mode, does not exist |j 
Therefore, demanding that the number of independent 
solutions (which is two) not reduce, one should start with 
p' 1 = 0{e~ 1 ). 



E. Slow-roll conditions and momentum constraint 

Interestingly, we can show that Bi vanishes when the 
slow-roll conditions arc satisfied. The slow-roll equations 
of motion for the canonical scalar fields are 



3M 2 H 2 ps V, 3ff V/v + Vi w 0. 



2, J 



(51) 



This implies that 



<Pn 



M 2 y, « (2e u - rnj)<p J N , (52) 



where the slow- roll parameters, e a b and r/at, are defined 
as 



M 2 P VjVj 

£ U = -TT^-T^-) VI J 



Then we find that the following relation is satisfied 
under the slow-roll condition, 

ipjfdiipx ps ip I N (2e IJ - r]u)d i ip J ps ip I NN d i ip I , (54) 

which yields Bi ps 0. In this sense, the slow-roll condi- 
tions are equivalent to the momentum constraint. 

What does this imply? This implies that the 6N for- 
malism happens to respect the momentum constraint, if 
the slow-roll conditions arc satisfied at the initial time, 

. This may provide a partial explanation as to why the 
6N formalism has been successful in computing £ for a 
wide variety of slow-roll inflation models. 



7 Assuming fl l = 0(e), Kodama and Hamazaki |1S(I also find this 
property from the momentum constraint. See their Eq. (5.13) 
and the argument given below it. 



V. CONCLUSION 

The necessary and sufficient condition for the validity 
of the SN formalism is that, with a suitable choice of 
coordinates, the perturbed Hamiltonian constraint and 
matter-field equations on large scales coincide with the 
corresponding unperturbed equations. 

That perturbed solutions in the long-wavelength limit 
can be obtained from unperturbed solutions was found 
and investigated by pioneering work in 1998 [T5L fl8l. . 
While their work was restricted to linear theory (and to 
quasi- linear theory [Hj]), we have extended their work 
to include non-linear (but still pcrturbative) perturba- 
tions. Such extension is also explored by la], who use the 
so-called "separate universe approach" [8j. As we have 
described in Sec. lIIIBl our starting point is more general 
than theirs. 

In this paper, using the flat gauge (ip = Cj = 0) and 
choosing the number of e-folds, N, as our time coordi- 
nates, we have shown that the perturbed Hamiltonian 
constraint and matter-field equations on large scales co- 
incide with the corresponding unperturbed equations, as 
long as (at least) the following conditions are satisfied: 

f. The unperturbed metric is given by a ho- 
mogeneous and isotropic Friedmann-Lcmaitre- 
Robertson- Walker metric. 

2. The final results for the curvature perturbation are 
not affected by decaying-mode terms, such as the 
shift vector, a time derivative of tensor perturba- 
tions to /a 2 (t), or the violation of the momentum 
constraint. 

3. Evolution of scalar-field perturbations outside the 
horizon can be treated using the lowest order or the 
next order of the gradient expansion. 

In order to show that the shift vector and a time deriva- 
tive of tensor perturbations are decaying modes, one 
needs two more conditions: 

4. Matter fields are given by scalar fields with an arbi- 
trary form of Lagrangian (whose anisotropic stress 
is of order 0(e 2 )) or, in the nth-order gradient ex- 
pansion where n = or 1, by some fluids with 
anisotropic stress of order 0(e n+1 ). 

5. Theory of gravitation determining the physics of 
inflation is Einstein's theory, or modified gravita- 
tional theories with can be transformed into Ein- 
stein's theory. 

The fourth and fifth conditions are sufficient conditions: 
we had to use Einstein's field equations to explicitly show 
that the fourth condition implies that the shift vector and 
a time derivative of tensor perturbations are decaying 
modes. It is possible that other theories of gravitation 
require different conditions for the shift vector and a time 
derivative of tensor perturbations to be decaying modes. 
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The discussion in th is p aper should also apply to 
vector-field models (see [20] for a review and references 
therein), as long as their anisotropic stress is of order 
0(e n+1 ) in the nth-order gradient expansion where n = 
or 1. 

The second condition is naturally expected in any in- 
flation scenarios. However, when the slow-roll conditions 
arc violated, it is known that decaying-modc solutions 
in the second order of the gradient expansion cannot 



be neglected in the power spectrum [21]. In such cases, 



since there are no gauges in which a Friedmann-Lemaitre- 
Robertson- Walker universe can be obtained in the second 
order of the gradient expansion, we can no longer use the 
5N formalism (22j. 

When all of the above conditions are satisfied, one can 
calculate £ using initial scalar-field perturbations com- 
puted in the flat gauge, and derivatives of the number 
of e-folds with respect to the initial values of the unper- 
turbed scalar fields, cpl, and their derivatives, yjy*. 



where n M = (l/a, — /3 l /aj is the unit vector normal to 
the ^-constant hypersurface, and V and D are the co- 
variant differential operators constructed by using g^ v 
and jij respectively. The dots denote time derivatives 
with respect to t. 

It is useful to decompose the extrinsic curvature Ki 3 
into a "trace" part K and a "trace free" part Aij as 



K ij = ^-K + a 2 e 2 ^A ij , 



(A4) 



where the indices of Ai 3 are raised/lowered by 7y, and 
A % i = is satisfied. 

By using the above notations, we write down Einstein's 
field equations. The Hamiltonian constraint is 



R (3) _ A l3 A^ + \k 2 = l^nV. 
The momentum constraint is 



(A5) 
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Appendix A: Einstein's field equations 

In this Appendix, we give Einstein's field equations in 
terms of the variables of the ADM formalism (see Eq. (1) 
for the ADM metric). We use Latin indices for the 3D 
spatial components running from 1 to 3, and Greek in- 
dices for the 4D spacetime components running from 
to 3. 

We decompose the three-space metric tensor as 

ll3 = a 2 (t)e 2 ^^ 3 , (Al) 

where a(t) is the scale factor. We define %j such that 
det[7ij] = 1; thus, det[7y] can be written as det^e' 1 ).^-] = 
e Tr I' 1 ] with Tr[h] = 0. We further decompose hij as 



(T) 

where h) 4 ' denotes a tensor mode and Ct has a scalar 
mode and a vector mode. 

The extrinsic curvature Ka is defined as 



d i c j + d j c i --s ij d k c k + h^\ 



(A2) 



Ki 



1 

2a 



(A3) 



(A6) 



(A7) 



The dynamical equation for ip is 

(d t - P%)iP + H = \(-aK + O^ 1 ). 
The dynamical equation for 7^- is 



(d t -p k d k )% j = ~2aA l3 +^ k d 1 l3 k +^ k d l (3 k ~-^d k l3\ 



which yields the dynamical equation for K: 

(d t - p k d k )K =a (.i,. ; .r + Ik 



k 

(A8) 



2M| 



(A9) 



as well as the dynamical equation for A^ : 



(d t ~ p k d k )A 
1 



R m _ 2p ( 3) \ _ ( D . D . a l^ DkD k, 



2 r 



+ a(KAij - 2A ik A K j ) + A lk d 3 f3 K + A jk d^ K - -A i3 d k p 



a 2 e W M 2 3 

(A10) 

Here, T,j — ^j-^ kl T k i is the anisotropic stress. 

The 3-dimensional Ricci scalar, RS 3 ', is constructed 
from jij, and M 2 is the reduced plank mass defined as 
l/87rG where G is the gravitational constant. 
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